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(N ; ABSTRACT 

' A sample of 123 radio sources that exhibit significant variations at 1.4 GHz on a seven year 

. base-hne has been created using FIRST VLA B-array data from 1995 and 2002 on a strip at 5 = 

' near the south Galactic cap. This sample spans the range of radio flux densities from 2 to 1000 

mjy. It presents both in size and radio flux density range a unique starting point for variability 
studies of galaxies and quasars harboring lower luminosity Active Galactic Nuclei (AGN). We 
find, by comparing our variable and non-variable control samples to the Sloan Digital Sky Survey 
the following: 1) The quasar fraction of both the variable and non- variable samples declines as 
^ , a function of declining radio flux density levels; 2) our variable sample contains a consistently 

higher fraction of quasars than the non-variable control sample, irrespective of radio flux; 3) the 
variable sources are almost twice as likely to be retrieved from the optical SDSS data than the 
' non-variable ones; 4) based on relative numbers, we estimate that quasars are about five times 

more likely to harbor a variable radio source than are galaxies; and 5) there does not appear to 
be any significant optical color offset between the two samples, even though the suggestive trend 
for sources to be bluer when variable has been detected before and may be real. This leads us to 
conclude that both radio variability and radio fiux density levels, in combination with accurate 
optical information, arc important discriminators in the study of (radio) variability of galaxies. 
The latter start to dominate the source counts below ^ 20mJy. In any case, variability appears 
to be an intrinsic property of radio sources, and is not limited to quasars. Radio variability at 
low flux density levels may offer a unique tool in AGN uniflcation studies. 

Subject headings: galaxies: active — galaxies: statistics — quasars: general 

1. Introduction vances in our understanding of radio source in- 

trinsic and extrinsic variability mechanisms have 
All is flux, nothing stays still; nothing endures ^een made by studying their amplitudes, charac- 

but changei. This seems to apply particularly teristic time-scales, and cross-correlations with ob- 

well to extragalactic radio sources. Important ad- servations at other wavelengths. Radio variability 

can be classified in three broad classes (Padrielli 



iReraclitus (540-480 BC) 



et al. 1987), based on their frequency depen- 
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dence: 1) variability only occurs at low frequen- 
cies (< 2GHz), 2) variability occurs both at low 
and high (> 5 GHz) frequencies, but is not corre- 
lated in time, and 3) low and high frequency vari- 
ability, but the variations are correlated (although 
they need not occur simultaneously). This classi- 
fication also conveniently separates the underlying 
variability mechanisms into "most likely extrinsic" 
(low frequency variations only), and "most likely 
intrinsic" (correlated broadband variability). 

The accepted model for the extrinsic variabil- 
ity of radio sources is refractive interstellar scin- 
tillation (ISS), which predominantly induccis varia- 
tions at low frequencies (< IGHz), with low ampli- 
tudes 2%), and with typical timescales of a few 
days (Blandford et al. 1986; Rickett 1986). Vari- 
ations of this type have been detected in multi- 
frequency studies (e.g., Mitchell et al. 1994; Lazio 
et al. 2001). Perhaps the strongest arguments 
for the external origin of these variations are cor- 
relations with Galactic latitude (e.g., Grcgorini 
et al. 1986; Gaensler et al. 2000), and an an- 
nual modulation due to earth's orbital motion 
(e.g., Lainela & Valtaoja 1993; Bondi et al. 1994; 
1996). Low- frequency variations with much longer 
timescales (~ 10 years), and/or larger amplitudes 
are thought to be source intrinsic, however (e.g., 
Rys & Machalski 1990 or, more recently, Gaensler 
& Hunstead 2000). 

Broadband radio variability, in which the vari- 
ations usually appear first and more strongly at 
higher frequencies, are thought to be intrinsic to 
the radio source. A model describing this behav- 
ior, in which shocks propagate along the radio jets, 
has been proposed by Marscher & Gear (1985). 
The model further implies that sources viewed 
close to the line of sight (i.e., flat-spectrum quasars 
and BL Lacs) should be more variable and have 
shorter variability timescales. For objects very 
close to the line of sight, relativistic beaming will 
become important, amplifying the variations and 
shortening the timescales to < 1 day (e.g.. Lister 
2001). Indeed, intraday variability has been read- 
ily detected in samples of compact, flat-spectrum 
sources. Some of the variability occurs on such 
short timescales that the apparent source bright- 
ness temperature (e.g., Gopal-Krishna et al. 1984) 
can be as high as T « lO^^K (e.g., Quirrenbach 
et al. 2000), well beyond the inverse Compton 
hmit of lO^^K (Kellermann & Pauhny-Toth 1969). 



Doppler boosting factors of up to a few hundred 
are needed to reconcile these numbers (the appar- 
ent T is proportional to the Doppler factor to the 
third power). These Doppler factors are actually 
at the high end of the observed range, and it has 
been suggested that these velocities are pattern 
speeds that do not reflect the Lorentz factors of 
the jet itself (e.g.. Lister et al. 2001). Nonetheless, 
intraday variability has been uniquely identified 
with the blazar class of objects (e.g., Quirrenbach 
et al. 1992; Lister et al. 2001). 

Comparatively little is known about (intrinsic) 
radio variability in non-blazar sources at low fiux 
density levels (~ 10 mJy). Most of the variability 
research has focused on either bright radio samples 
(e.g., Grcgorini et al. 1986, S > 0.4 Jy ; Lister et 
al. 2001, S > 0.4 Jy; Aller et al. 2003, S > 1.3 
Jy), or small samples at much lower flux levels 
(a few lO's of mJy - e.g., Rys & Machalski 1990; 
Bondi et al. 1996; Riley et al. 1998). Recently, 
Carilli et al. (2003) used deep multi-epoch VLA^ 
observations of the Lockman Hole area to study 
variability at the 0.1 mJy level, resulting in just 9 
variable sources. 

In this paper, we will use two epoch observa- 
tions of the FIRST Zero-Dec strip (Sect. 2) of 
about 9000 radio sources to expand significantly 
the number of variable sources down to flux den- 
sity levels of '--^ 1 mJy. Unlike much brighter ra- 
dio samples, the typical radio source population 
at these levels is dominated by star-forming galax- 
ies and steep-spectrum systems, with only a small 
percentage of flat-spectrum AGN (e.g., Windhorst 
et al. 1999; Richards et al. 1999; de Vries et al. 
2002). With our sample, therefore, we are in an 
excellent position to investigate the radio variabil- 
ity properties of lower radio luminosity galaxies. 

Interstellar scintillation is not expected to 
contribute signiflcantly to the variability in our 
sources, since our observations were carried out 
at 1.4 GHz, near the minimum in expected ISS 
(e.g., Padrielh et al. 1987; Mitchell et al. 1994). 
Furthermore, given our variability criterion (see 
Sect. 3), only sources that vary by significant 
amounts, especially toward the lower flux levels, 

^The Very Large Array is a facility of the National Radio 
Astronomy Observatory (NRAO) , which is operated by the 
Associated Universities, Inc., under cooperative agreement 
with the National Science Foundation. 
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Fig. 1. — FIRST Zero-Dec survey area. Note the non-uniform source density cfose to the survey boundaries 
(due to the very different axis scales this is not obvious for the RA edges). 



are included. These variation levels exceed the 
typical ISS variations of a few percent by a large 
factor. 

The paper is outlined as follows. In the next 
two sections, the radio observations and sample 
selection are described. In Sect. 5 we discuss the 
optical properties of the variable sample, and com- 
pare these against a non- variable control sample. 

2. Observations 

Our radio imaging data have been taken as 
part of the Faint Images of the Radio Sky at 
Twenty centimeters survey (FIRST, see Becker et 
al. 1995 for a detailed description) conducted with 
the VLA in B-array between 1993 and 2002. The 
observing strategy of the FIRST survey provides 
a limited capability to search for variability, al- 
though since adjacent fields have significant over- 
lap, most sources were observed more than once. 
For the most part, data were collected in strips 
of constant declination so adjacent fields in the 
east-west direction are sensitive to variations on 
a time scale of 3 minutes. FIRST observing runs 
were usually separated by multiples of 24 hours 



so adjacent fields in the north-south direction are 
sensitive to variations on timescales of 1 to 5 days. 
This led, for instance, to the discovery of a number 
of radio- variable stars (Helfand et al. 1999). 

A search for variability between the FIRST 
survey and the NRAO VLA Sky Survey (NVSS, 
Condon et al. 1998) is also possible. This is, 
however, complicated by the difference in an- 
gular resolution: the FIRST survey can resolve 
out flux seen by the NVSS. This results in so 
many NVSS sources appearing brighter than their 
FIRST counterparts that it is necessary to restrict 
the analysis to sources which appear brighter in 
FIRST (see Sect. 3.2), making a straightforward 
flux comparison to test for variability difficult. 

The best opportunity to search for long term 
variability has been afford by our decision to re- 
observe the south Galactic cap (SGC) FIRST Zero 
Dec strip (hereafter FZDB) during the summer of 
2002 (the initial set of observation were made in 
1995, which wiU be referred to as FZDA). The 
observations were made for the dual purpose of a 
quality control test of the FIRST survey, as well as 
a search for variability. The zero-dec strip was cho- 



3 



T-i-rj 




J ' ' ' ' ' 'I 

1 10 

Radio Flux Density [mJy] 

Fig. 2.— Fractions of the FZDA and FZDB cat- 
alogs (green and red lines respectively) that are 
present in the combined catalog. The histograms 
are for point sources only, and the matching has 
been done to within 3". Note the high fractions for 
sources with flux densities exceeding 2 mJy (98.0% 
on average) . It is also apparent that some fraction 
of FZDB sources below 2 mJy are spurious, more 
so than for the FZDA catalog. Each bin represents 
500 sources, so the closer spacing toward lower 
flux density levels reflects the increasing number 
of faint objects. 

sen for this purpose because it has been observed 
by the Sloan Digital Sky Survey as well (SDSS, see 
Stoughton et al. 2002 for the Early Data Release, 
and Abazajian et al. 2003 for Data Release 1 up- 
dates). A total of 720 fields were re-observed cov- 
ering 120.2 square degrees and encompassing 9086 
radio objects in the FZDB catalog. The FZDB 
survey area is shown in Fig. 1. 

2.1. Matching FZDB with FIRST 

The FIRST catalog has a uniform sensitivity 
across the surveyed zero-dec strip area, whereas 
the FZDB has not (Fig. 1). This makes a di- 
rect comparison of their relative completeness lim- 
its and source surface densities difficult. With an 
average surface density of ~ 90 radio sources per 
square degree White et al. (1997), the FIRST cata- 



log should contain about 10 800 objects within the 
nominal zero-dec strip survey area. However, this 
number cannot be directly compared to the FZDB 
count of 9086, due to the latter's non-uniform sur- 
vey coverage. The best way to compare the two 
catalogs is to create an FZDA catalog (derived 
from FIRST) which has its source surface den- 
sity tapering off toward the survey edges, just like 
the FZDB catalog. The simplest way of obtaining 
such a catalog is to list the unique sources con- 
tained in FIRST around the FZDB catalog posi- 
tions. This makes sense, since we are interested in 
objects present in both catalogs in the first place. 
The total number of sources in such an FZDA cat- 
alog vary with the radius up to which the FZDB 
positions are matched: 6605 sources within 3", 
7343 within 1', and 11778 within 10'. We will 
use the first two lists to compare the relative com- 
pleteness of the FZDB catalog. The 6605 entry 
list represents our matched catalog from which we 
will select our variable sample (Sect. 3.1), and the 
7343 constitutes our comparison FZDA sample. 

The combined catalog uses a matching radius 
of 3", which is about half of the FWHM of the 
B-array beam. This automatically excludes the 
fraction of extended sources for which the source 
detection program assigned different components 
and/or different positions though. As a conse- 
quence, the number of 6605 is considerably lower 
than the original FIRST surface density (even 
lower than expected based on the non-uniform 
coverage). To be consistent, we exclude all ex- 
tended components^ from our combined sample. 
There are good scientific reasons for this as well. 
Extended radio emission is most likely due to 
lobe (radio hotspot) emission which cannot vary 
on the short timescales we are sampling here 
(7/(1-1-2:) years). Since we are interested in shorter 
timescales, we will have to study those sources in 
which the (feeding) processes close to the AGN 
itself vary. These sources are consequently very 
small, and unresolved by our imaging data (with 
its resolution of 5 arcseconds). The resulting 

^We consider a source / component extended if it has a 
peak flux less than 70% of its integrated flux. It should be 
noted that we are referring to components, so, for instance, 
it is possible that a single extended radio source can be 
resolved into a number of unresolved components, each of 
which would end up individually in our combined sample; 
see Sect 4. 
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reduced sample sizes have been listed in Table 1. 



2.5 



2.1.1. Sources unique to one catalog 
phenomena 



Transient 



Figure 2 illustrates how well the FZDA and 
FZDB compare. At the higher flux density levels 
we expect to find almost all of the point sources in 
the combined catalog. This is the case for sources 
with flux densities exceeding 2 mJy. For both 
the FZDA and FZDB catalogs, over 98% of those 
sources end up in the combined catalog. None of 
the remaining catalog-unique entries turned out 
to be real on close inspection: they were either 
resolved sources for which a particular component 
was assigned a differing position between FZDA 
and FDZB (separated by more than 3"), or the 
catalog entry turned out to be an imaging artifact 
that made it into one but not both catalogs (e.g., 
ringing around very bright objects). 

The situation is quite different for fainter 
sources, however. An increasing fraction of faint 
objects are only present in either catalog, which 
indicates a systematic effect rather than genuine 
source variability. Disentangling the truly vari- 
able sources from the random noise at these flux 
density levels can only be achieved by obtain- 
ing deeper surveys. Indeed, none of these faint 
sources end up meeting our variability criterion 
(Sect. 3.1), which is based on the survey sensitiv- 
ity level. All the variable sources that do meet the 
criterion have flux densities > 2mJy (see Fig. 5), 
and are in the well-matched part of Fig. 2. 

So, in summary, transient phenomena (e.g.. 
Gamma Ray Burst afterglows) have not been 
found for flux density limits exceeding 2 mJy. 
Fainter events may be present, but we do not 
have any way of assessing their reality. 

3. Variable sample selection 

Since we are interested in radio source vari- 
ability, we begin by deflning our Variability Ratio 
(hereafter VR) as: 



VR(i) 



(1) 



for the i^^ source in our catalog. The indices A 
and B denote the measurement epochs 1995 and 
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Fig. 3. — Relative flux density ratios for our point 
source sample as function of S/N. The flux density 
ratio is deflned here as larger flux density divided 
by smaller flux, irrespective of the epoch (hence 
it is always > 1). The red points are for the inte- 
grated flux density values, and the green points for 
just the fitted peak values. As can be seen by the 
overplotted medians (purple for integrated, blue 
for peak fiux densities), the former has a larger 
scatter, especially toward lower flux density lev- 
els. 



2002. Because the source brightness does not en- 
ter into this equation, we cannot use this ratio 
at face value: toward the lower flux density lev- 
els the (sky background) noise will account for 
a significant fraction of the "variability" in such 
sources. We therefore have to weight the VR value 
of each source individually based on its Signal-to- 
Noise Ratio (S /N) and local noise values (see next 
section). 

Since we are only considering unresolved sources, 
ideally their (fitted) peak fiux density and inte- 
grated flux density should be the same. This is, 
however, not the case. In fact, the integrated flux 
density values are found to have larger intrinsic 
scatter than the fltted peak flux values. Both 
the peak and integrated flux density are based on 
Gaussian flts to the actual data, but unlike the 
peak fit, the integrated value is also based on flts 
to the major and minor axes of the object. Noise 
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Fig. 4. — Flux density variations for unresolved 
sources between the two epoch as function of 
signal-to-noise. The peak flux density value of the 
most recent epoch has been divided by the 1995 
data. The yellow lines represent the 4ct variation 
limits, as given by Eqn. 4. The highlighted green 
crosses are sources that meet our variability crite- 
rion, the blue crosses are considered non-variable. 
The solid green line marks the local median vari- 
ation value, and serves as an indication of the rel- 
ative calibration between the two epochs. This 
line has been replotted for better clarity on a lin- 
ear scale (indicated on the right vertical axis) as a 
solid red hue. 

in the measurement of the axes therefore increases 
the noise in the integrated flux density compared 
to the peak value. In Fig. 3 we have plotted the 
relative flux density difference ratio (defined as 
the larger flux density divided by the smaller flux 
density) for both the source integrated and peak 
flux density values. 

It is clear from this figure that the scatter in 
the peak flux measurements (the solid blue me- 
dian line) is considerably less than that of the in- 
tegrated flux density ratios (as indicated by the 
solid purple median line). To minimize the num- 
ber of variable false positives, we only consider the 
source peak flux density values. 
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Fig. 5. — The same plot as Fig. 4, but with the 
adjusted FZDB flux densities: 1.0116 x {Sb + 
0.090 mJy). This straightens out the median vari- 
ability curve and evens out the number of positive 
and negative variable sources. 

3.1. Variability criterion 

Following Rengelink et al. (1997), the relative 
flux density errors can be written as: 



(2) 



This equation reflects the two components of the 
measurement error, with C'l due to a constant sys- 
tematic error, and C2 dependent on the S/N ratio. 
This expression can be used to estimate the sig- 
niflcance of radio source variation, as defined by 
the quotient of the new and old flux densities. If 
one defines the flux density and noise values in the 
first epoch as Sa^cta, and in the second as 83,0-3, 
the 4cr variability envelopes are then given by: 



VR4 



S3 {i±^^ci + ci{^y' 

SA{lT^^Cl+Cl{^^f 



(3) 



which can be further reduced by substituting Sa = 
Sb = Sf^y-^Q, (Ja = 0-3 = 0.15 mJy, Ci = 0.01, and 
C2 — 1.3. The a values have been measured to be 
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on average 0.150 and 0.148 mJy for the 1995 (A) 
and 2002 (B) epochs, respectively. The Ca value 
has been chosen to match the value of Rengelink 
et al. (1997) and dc Vrics ct al. (2002), and seems 
fairly robust across different instruments and sur- 
vey setups. The C\ = 0.01 value (which is lower 
than the Ci = 0.04 of the latter two papers), is 
set by the high S/N end of the source distribu- 
tion. Its adopted value is a good match between 
including / excluding high S/N variable sources: 
the fraction of variable sources as function of S/N 
is kept acceptably flat (see Figs. 4 and 5, with S/N 
defined as (Sa + Sb)/{(ta + fe))- 

The upper and lower Aa envelopes can now be 
expressed (as a function of S/N) as: 



yj^upper 



4x 



4x 



yj^k>wer 



= (VR4" 



with: X = VO.OI^ 



/ 1.3 



VS/N 



(4) 



Equation 4 forms the basis of our variability 
criterion. Any source having a VR exceeding ei- 
ther 4(7 envelope is considered genuinely variable. 
This yields 146 variable sources, of which 58 are 
brighter in the 2002 epoch, and 88 are fainter 
(Fig. 4). 

The solid red curve in Fig. 4 represents the me- 
dian flux density ratio as function of S/N, plotted 
on a linear scale to improve clarity (as opposed to 
the logarithmic equivalent in green). As can be 
gathered from this curve, the medians quickly fall 
below 1 as the S/N decreases, with a maximum 
deviation of 8% at S/N ~ 10. Most of the offsets 
are within ~ 2%, however. Even though the off- 
sets are small, there is a systematic trend which 
is affecting our variability sample selection, some- 
thing we wish to avoid. If uncorrected, it implies 
that fainter sources are more likely to fade between 
1995 and 2002 than brighter sources, something 
which would be hard to explain in any variability 
scenario. We have two parameters to play with: 
a small (additive) zero-point correction which will 
straighten out the median curve, and a multiplica- 
tive sensitivity correction which will raise or lower 
the curve along the y-axis. The best correction 
(with the resulting median curve having the small- 
est variations along the VRc = 1 line) is obtained 



by adjusting the FZDB flux densities by: 



_ 1.0116 X (g(i)B + 0.090 mJy) 
VRc(^) - (5) 

which amounts to a 90 /iJy zero-point offset, and 
a 1.16% sensitivity correction. The 90 /iJy offset 
is small compared to both the 150 /iJy rms noise 
value and the 250 /zJy CLEAN bias (a flux den- 
sity correction applied per beam for the FIRST 
survey, see White et al. 1997). And the 1.16% 
factor is well within the nominal ~ 5% systematic 
uncertainty in the flux density scale. 

After applying these corrections (Fig. 5), we 
end up with 128 candidate variable sources / com- 
ponents (corresponding to 2.0% of the sample), 
listed in Table 2, of which 70 are brighter in the 
2002 epoch, and 58 fainter. 

3.2. NVSS comparison 

As a reality check, we compared the FZDB 

catalog against the NVSS catalog to see whether 
the same sources are deemed variable, using ba- 
sically an identical variability criterion. There 
are, however, a few differences. The NVSS sur- 
vey has a typical background noise of 0.45 mJy 
(Condon et al. 1998), instead of the 0.15 mJy 
for the FIRST survey. Another complicating fac- 
tor is the resolution difference between the NVSS 
and FZDB. Sources that are not resolved by the 
45" NVSS beam are often resolved by the FZDB, 
with each component only a fraction of the NVSS 
total flux density. We tried to minimize this ef- 
fect by first excluding resolved NVSS sources out- 
right from the correlation, and secondly using the 
FZDB integrated (instead of peak) flux density 
measurements. Furthermore, we used a very re- 
strictive 5" matching radius in another attempt to 
just match unresolved sources both in the FZDB 
and NVSS. The resulting matched sample contains 
2201 sources. 

Of these 2201 sources, 10 are brighter at the 
FZDB epoch, and 91 are fainter. This is after ap- 
plying a similar FZDB flux density correction as in 
the previous section (albeit with slightly different 
constants) . Clearly, there is still some flux density 
being resolved out by the longer base-line FZDB 
observations compared to the NVSS, which results 
in quite a few bogus variable sources. 
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This implies wc cannot really compare the 
FZDB-NVSS faint variable sources to Table 2, 
since an unknown fraction of them will be spuri- 
ous. The 10 FZDB-NVSS bright variable sources 
on the other hand, are not be affected by this. In- 
deed, 9 out of the 10 sources are recovered in Ta- 
ble 2 (separately footnoted) , with the sole "newly" 
identified variable source FIRST J222646+005210. 
Upon closer inspection, this source appears to be 
stable between FZDA and FZDB (653 ± 1 mJy for 
both epochs). Its NVSS flux density of 615 mJy is 
just within the quoted ~ 5% systematic flux den- 
sity uncertainty between the two surveys. Based 
on this, we therefore conclude that our variabil- 
ity criterion is robust, and results in a sample of 
truly variable sources. This does tacitly assume 
that the cataloged radio flux densities are always 
correct. However, as the next section will show, 
this is actually not valid for some cases. 

4. Radio morphology 

Our variable sample actually contains 128 vari- 
able, unresolved, components. While most of these 
component entries turn out to be single (unre- 
solved) radio pointsources, some of them are part 
of a multi-component radio source. The varying 
part is then presumably associated with the AGN 
/ radio-core position. However, there may still be 
some cases where the FZDA and FZDB cataloging 
software divided up radio flux differently between 
adjacent components, introducing artificial vari- 
ability. It is therefore important to check all the 
multi-component radio sources in our sample (28) 
individually. 

In order to characterize the overall radio prop- 
erties of our sample, we start by defining 6 mor- 
phological categories: 1) Single pointsource (PS 
in Column 11 of Table 1); 2) Pointsource with an 
elongation, or linear feature (CJ, "core-jet"); 3) 
Variable core with two non- variable lobes (CL); 4) 
Core embedded in diffuse halo (CH); 5) Complex, 
multi-component structure (CX); and 6) Possible 
hotspot variability (HS). 

Each of the sources in Table 1 have been cate- 
gorized accordingly. The relative number break- 
down of the sample is: 1) 100 (78.1%); 2) 14 
(10.9%); 3) 9 (7.0%); 4) 1 (0.8%); 5) 1 (0.8%); 
and 6) 3 (2.3%). Except for the last two classes 
(containing 4 sources, or 3.1% of the total sample). 



all of the variability is consistent with it being due 
to the AGN / core component. We will discuss a 
few selected sources in more detail below. Radio 
maps of all the extended variable radio sources can 
be found in Fig. 10. 

J000257-002447 (CH) - This source contains 
a variable pointsource component embedded in an 
extended radio halo. The radio surface brightness 
of this halo remains constant between the FZDA 
and FZDB epochs, whereas the peak flux den- 
sity of the central pointsource varies by more than 
10%. 

J001800-K000313 (HS) - The North-Eastern 
lobe of this double radio source is slightly resolved, 
just enough so that it did not make it into our vari- 
able component list. As it turns out, both lobes 
/ hotspots appear to be variable. The total flux 
density for the radio source increased by 13.6% 
between FZDA and FZDB, from 62.14 to 70.56 
mJy. Unfortunately, two unrelated sources close 
to this source also "varied" by ^ 10%, casting se- 
rious doubt on its reality. The source is very close 
to the edge of its particular survey fleld, which may 
be the reason for the flux density discrepancy. It 
is taken out of further consideration, and has been 
placed separately in Table 1. 
J003246-001917 (HS) Upon close inspection 
of both the FZDA and FZDB catalogs, we have 
to come to the conclusion that this source is not 
variable. The problems stems from the fact that 
FZDA lists 5 components for this source, and 
FZDB only 4. This leads to quite different flux 
density assignments to the North- Western lobe. 
The source is taken out of further consideration. 
J012213-001801 (CL) - Both the Eastern and 
Western lobes are constant to within ~ 1% be- 
tween FZDA and FZDB. Flux density values for 
the bright Eastern lobe are 123.84 and 122.61 mJy, 
a variation of 1.0%. The core, on the other hand, 
brightened over the same time period from 348.28 
to 404.79 mJy. 

J020234+000301 (CL) - The lobes do not vary 

beyond the 1% level. The core, however, fades 
from its FZDA flux of 44.01 mJy to 35.92 (an 
18.4% decrease). 

J021840-001515 (CX) - This is the only source 
in our sample with a hard to categorize radio 
morphology. The marked component position 
(Fig. 10) coincides with a, z = 1.171 quasar. The 
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two components immediately to the North do not 
have optical counterparts, though there appears 
to be a faint object in the SDSS image that is 
situated in between the lobes. All of the radio 
components fade between FZDA and FZDB by on 
average 16%. A nearby, unrelated, radio source 
observed in the same field remains constant, how- 
ever. The component associated with the quasar 
varies the most: from 14.12 mJy down to 10.92 
mJy (a decline of 22.7%). At the moment it is not 
clear to us how the various components relate to 
each other, if at all. 

J021202-002750 (CL) The core component 
faded by at least 4 mJy (~ 8%) relative to the 
lobes. 

J212058+000612 (HS) - As is the case for 

J003246-001917, this source is not variable. The 
South-Eastern lobe is constant, and the perceived 
variability in the North- Western lobe is due to dif- 
fering FZDA and FZDB component flux density 
assignments. This source is also taken out of con- 
sideration. 

J212955+000758 (CL) - The "variability" is 
due to the FZDA having 3 components, whereas 
FZDB lists 4. The summed fluxes do not vary. 
J220017-000133 (CL) - The core component in 
this source brightened (by 20%), so we are inclined 
to regard this source as variable even though the 
number of fitted components is different between 
the FZDA and FZDB. This faint extra component 
in the FZDA catalog (just to the North of the core) 
docs not appear to contribute anything to the total 
flux of the core in the FZDB catalog (where this 
component is absent), based on the fitting param- 
eters. 

J222729-h000522 (CL) - The core variability is 
17.2% (relative to the lobes), with the flux declin- 
ing from 97.51 mJy to 80.69 mJy. 
J235050-002848 (CL The core flux increased, 
whereas the lobes remained more or less constant. 
J235828+003934 (CL - All of the components 
faded by about the same amount, casting serious 
doubt on the reality of core variability. We have 
taken this source out of further consideration. 

So, after close inspection, we found all 3 poten- 
tial hotspot variable sources (HS) and 2 core-lobe 
(CL) sources to be non-variable. They have been 
excluded from our subsequent studies, but are left 
in Table 1 and Fig. 10 as reminders of how careful 



one has to be using automated selection criteria. 

5. Optical identifications - Detection rates 

We next investigate the possible connection be- 
tween optical properties and the presence of ra- 
dio variability. Variability correlations in the op- 
tical and radio have been proposed before (e.g., 
Tornikoski et al. 1994; Hanski et al. 2002), with 
the highly beamed blazars as the clearest examples 
of sources which exhibit variability across the elec- 
tromagnetic spectrum (e.g., Bregman et al. 1990). 
The first step would be to correlate our variable 
sample with known optical catalogs, and compare 
the results against a sample of non-variable ra- 
dio sources. We therefore created 12 sets of 123 
non-variable sources (i.e., each control sample has 
the same number of elements as the variable sam- 
ple). To ensure non- variability, we only selected 
sources inside la variation curves (as opposed to 
sources outside 4cr variations. Fig. 5), while match- 
ing the radio flux density distribution function of 
the variable sample. This was done by matching 
the number of non-variable sources to the amomit 
of variable sources per S/N bin (binned logarith- 
mically from 1.0 to 3.4, with binsize 0.4 dex). 
Given that there are only 1922 sources (out of 
our initial sample of 6605) that are less variable 
than these la limits, 12 sets of control samples are 
about as many as one can create without starting 
to have a signiflcant fraction of sources shared be- 
tween them. 

All of our FZDA/FZDB survey area is covered 
by the Automatic Plate Measuring (APM) ma- 
chine catalog. The APM facility (in Cambridge, 
UK) lists identiflcations and positions based on 
scanned UK and POSS II Schmidt plates, cover- 
ing currently more than 15, 000 deg^ of sky. 

5.1. APM matches 

Our sample of 123 variable sources has been 
matched against the APM catalogs with a 3" 
matching radius. While this does discriminate 
against possible off-nuclear radio sources, it is nec- 
essary to limit the number of chance matches, 
given the density of objects in the APM catalog. 
The matching results are listed in Table 3. What 
is immediately striking is the difference in detec- 
tion rates between the variable and non-variable 
samples. About half of the variable sources are 
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10 15 SO 

(B+R)/2 [mag] 

Fig. 6. — Cumulative magnitude distribution of 
radio sources which are identified by the APM (66 
variable sources in total). The solid green line 
is for our variable sample, the dashed red line 
represents the distribution for the combined 12 
non-variable (less than la) control samples (262 
sources in total) . Both distributions are identical, 
based on a KS-test (using 0.25 magnitude bins, for 
a total of 56). 

retrieved from the APM catalogs, but only 1 in 
6 of the non-variable sources have optical coun- 
terparts. The rms variation between the 12 con- 
trol samples is 2.8%, resulting in a 13a detec- 
tion rate difference. Clearly, variable sources are 
easier to detect not only because on average they 
are brighter, but also due to a higher fraction of 
quasars^ (see Sect. 5.2). The brightness difference 
is not large; actually a cumulative distribution plot 
against mean magnitude (Fig. 6) shows how com- 
parable their distributions are. Based on a KS- 
test, these two distributions do not differ to any 
level of significance. 

Fortunately, a significant fraction ('^ 60%) of 
the survey is covered by the SDSS Data Release 
1 (DRl), so wc can use these deeper and more 
uniform imaging data to investigate this further. 

''without spectroscopic confirmation, most of tfiese objects 
sfiould be considered quasar candidates. 



5.2. SDSS 

The SDSS data also show a higher fraction of 
optical identifications for variable sources. Since 
the latter data are deeper than the photographic 
plates the APM catalogs are based on, we identify 
82% of the variable sources, and on average 38% of 
the non- variable ones (Table 4). These values are 
up from the APM values of 54% and 18%, respec- 
tively. As mentioned earlier, one way of explain- 
ing the higher match rate for the variable sample 
is that it may contain a higher fraction of quasars. 
Since these objects are typically unresolved, they 
are easier to detect at any given magnitude than 
resolved galaxies (of the same integrated magni- 
tude). Table 4 lists the morphological breakdown 
of each sample, based on the SDSS classification. 
With a relative stellar fraction of 53%, the variable 
sample contains significantly more quasar candi- 
dates than the non- variable control sample (on av- 
erage 19%). 

The results of Table 4 can also be used to cal- 
culate how much more quasars are likely to be 
variable compared to galaxies. This is a function 
of radio flux density (Fig. 7), but the low num- 
ber of sources (per bin) only allow for an overall 
ratio. We assume the non-variable sample num- 
bers (19% of the radio population are quasars, 
81% are galaxies) are representative of the gen- 
eral population. These numbers compare very well 
to the ones quoted by Ivezic et al. (2002) for a 
much larger sample of SDSS sources detected with 
FIRST (83% of those are galaxies). Now, if we in- 
troduce the terms quasar variability rate QVR and 
galaxy variability rate GVR (both defined as the 
fraction of quasars / galaxies that are variable), 
and an initial sample size of P, we get: 

0.19 X P X QVR = # of quasars in var. sample 
0.81 X P X GVR = # of galaxies in var. sample (6) 

This means that the variable sample will contain 
in total (0.19 xPxQVR-|-0.81xPx GVR) variable 
sources. We know the relative quasar and galaxy 
fractions to be 53% and 47%. This leads to (after 
dividing out P): 

0.19 X QVR = 
0.53 X (0.19 X QVR + 0.81 x GVR) ^ 
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Fig. 7. — Quasar population fraction as a func- 
tion of radio flux density. The solid green line 
denotes the quasar fraction for our sample (64 ob- 
jects in total, 16 per bin), and the dashed red hue 
indicates the fraction for our non-variable control 
sample (341 sources in total, 60 per bin). The 
error-bars are given by the formal binomial distri- 
bution error: a{f) = (/(I - /)/A^)°^ 



_ 0.53 X 0.81 
^^^-0.47x0.19 

= 4.8 GVR 



(7) 



If one assumes a ^ 6% uncertainty in the rela- 
tive fractions (Table 4), this relative value ranges 
between 3.4 and 7.5. In other words, not surpris- 
ingly, quasars are on average 5 times as variable 
(in the radio) than are galaxies. 

5.3. Classification as function of radio flux 
density 

The radio source population changes as a func- 
tion of radio flux density (e.g., Windhorst et al. 
1999; Richards et al. 1999; de Vries et al. 2002), 
from mainly AGN-dominated sources at higher 
flux density levels down to starburst dominated 
sources at the sub-mJy level. This will affect our 
variable sample as well, and we expect to see a de- 
cline of the quasar fraction as function of lower flux 
density limits. This is nicely illustrated in Fig. 7, 
which plots the relative quasar fraction with flux 



density. The solid green line represents our vari- 
able sample, and the dashed red line is the mean of 
our 12 control samples. The bins have a constant 
number of members (except for the last bin in the 
control sample), hence the variable spacing. The 
binning of the control samples has been adjusted 
to approximately match the binning of the actual 
sample (4 bins of 16 each for the variable sample, 
and 5 bins of 60 plus 1 bin of 41 members for the 
control samples). 

While it is not known how much uncertainty 
is present in the SDSS classification (and hence 
quasar fraction) toward the faintest magnitudes 
(see Sect. 5.3.1 below), it is apparent from Fig. 7 
that the quasar fractions for the variable and non- 
variable samples differ significantly (the formal 
fraction errors are 5% for the control sample, and 
12% for the variable sample). This leads us to two 
robust conclusions based on this figure: 1) a radio 
variable sample will contain more quasars than a 
non- variable sample, irrespective of radio flux den- 
sity, and 2) the relative contribution of quasars to 
the sample declines as a function of declining flux 
density: below ^ 20mJy, galaxies account for the 
majority of the variable sources. 

5.3.1. Classification reliability 

The last section critically hinges on the reli- 
ability of the SDSS galaxy versus quasar classi- 
fication scheme. If for one reason, most of the 
excess quasar classifications in our variable sam- 
ple were made toward the very faint end of the 
magnitude range, one might become suspicious 
about the reality of the effect. A direct test is 
plotting the quasar fraction as function of optical 
(R-band in this case) magnitude, instead of radio 
flux density as in Fig. 7. This has been done in 
Fig. 8, with again our variable sample in green 
and the control sample in red. It is reassuring to 
see that most of the excess quasar identifications 
in the variable sample are made for magnitudes 
brighter than R ~ 21. Stoughton et al. (2002) 
estimate a 95% confidence level for classifications 
up to this magnitude. It becomes much less cer- 
tain for fainter sources. Without doubt some of 
the current classifications are wrong, but there is 
no indication that our consistently higher quasar 
fraction in the variable sample, both as function 
of radio flux density and optical magnitude, is sig- 
nificantly affected by this. 
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Fig. 9. — SDSS (r— z) versus (u— r) color plots. The variable sample is represented in green, with quasars 
and galaxies in the left and right panels respectively. The control sample is presented in red. The crosses 
mark the sample mean colors and their la errors along each axis (values are given in Table 5). The optical 
color offset between the variable and non-variable samples are significant at the ^ 2a level. This trend for 
variable sources to become bluer during an outburst has been seen before (e.g., Trevese & Vagnetti 2002). 
Note also that the the quasars are generally bluer (i.e., lower u— r and r— z colors) than the galaxies. 



5.4. Optical colors 

We have established that there is a significant 
difference between the variable and non-variable 
radio samples in terms of their optical counter- 
parts. The former have a higher matching rate 
with optical identifications (given a fixed detec- 
tion limit), which is at least in part ascribed to 
the higher fraction of quasars. We have already 
seen that the magnitude distributions are not sig- 
nificantly different (Fig. 6), however, this does not 
rule out that significant optical color differences 
may be present between the samples. The clear- 
est difference between galaxies and quasars (at 
least in the restframe) is provided by the (u— r) 
vs (r— z) color plot (see Fig 13 in Stoughton et al. 
2002 for intermediate color-color plots). This al- 
lows for sampling of the 4000A break in galaxies 
up to a redshift of about unity. This break will 
offset the galaxies from the quasars, which usually 
exhibit fairly flat, powerlaw continua (see Fig. 3 
in de Vries et al. 2003). Our data are presented 
in Fig. 9. The sample means and their la errors 
(indicated by the crosses in Fig. 9) are presented 
in Table 5. Both the variable quasar and galaxy 
distributions in the color-color plane are (statisti- 



cally) distinct from the non-variable control sam- 
ple. This trend for both the galaxies and quasars 
to become bluer when they vary has been seen be- 
fore (e.g., Giveon et al. 1999; Trevese & Vagnetti 
2002; Vagnetti et al. 2003; de Vries et al. 2003). 
Still, the color spread among the sources in the 
sample is much larger than the individual color 
change as a source goes through an outburst. In- 
deed, the studies of Giveon, Trevese, Vagnetti and 
their respective collaborators all rely on following 
individual objects through their outbursts, instead 
of using a statistical sample. 

The presence of the mean optical color offsets 
in our radio variable sample implies two things: 
1) on average, radio variable sources also exhibit 
optical (color) variations, and 2) the optical vari- 
ation time-scales are shorter than our 7 year base- 
line, since the SDSS data were not taken at the 
epoch of radio outburst. However, optical colors 
for an individual source, while useful in differen- 
tiating between quasars and galaxies, do not pro- 
vide an effective selection mechanism for (radio) 
variability. Nonetheless, as we have seen in the 
previous sections, the presence of radio variability 
has a clear impact on the sample's morphological 
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18 20 
SDSS R-band [mag] 

Fig. 8. — Quasar population fraction as function 
of R-band magnitude. Like in Fig. 7, the solid 
green line denotes the quasar fraction for our vari- 
able sample, and the dashed red line is for the 
non- variable control sample. We applied slightly 
larger bins here (20 and 90 objects for the vari- 
able and control samples). Note that most of the 
excess quasars in the variable sample are toward 
the brighter magnitudes (R<21), and can be con- 
sidered reliably classified. Error-bars as defined in 
Fig. 7. 



makeup. 

6. Summary 

We have created a sample of 123 variable radio 
sources using two epoch observations of a zero-dec 
strip toward the south Galactic cap. This sam- 
ple spans the range of radio flux densities from 
~ 2 to 1000 mJy. It presents both in size and ra- 
dio flux density coverage a unique starting point 
for variability studies of more normal, less AGN- 
dominated galaxies, especially toward the lower 
flux density limits. We compared both our vari- 
able and non- variable samples to the Sloan Digital 
Sky Survey optical data. 

We found that the quasar fraction of the sam- 
ple sharply declines as a function of declining radio 
flux density levels. This is consistent with earlier 
findings that the radio source population demo- 



graphics change as one samples at progressively 
lower flux density levels: AGN-dominated systems 
tend to be found at the brighter radio flux density 
levels (> 10 mJy), whereas star-forming and nor- 
mal galaxies dominate the counts at sub-mJy flux 
density levels (at 1.4 GHz, e.g., Windhorst 2003). 
Our variable sample contains a consistently higher 
fraction of quasars than the non-variable control 
sample, independent of radio flux density. While 
this explains part of our almost 2 x higher optical 
matching rate of the variable sample compared to 
the non-variable one (quasars are easier to detect 
at a given brightness limit than galaxies), it does 
imply that our variable sample contains on aver- 
age slightly brighter sources (though not signifi- 
cantly so, see Sect. 5.1). Based on relative num- 
ber statistics, we estimate that quasars are about 5 
times more likely to harbor a variable radio source 
than galaxies. However, at flux density levels < 20 
mJy, the majority of (radio) variable sources are 
identified as galaxies. And finally, galaxies and 
quasars that harbor a variable radio source ex- 
hibit, on average, bluer optical colors than hosts 
of non- variable sources. 

All of this underlines the fact that both galaxies 
and quasars can harbor variable radio sources, al- 
beit at different occurrence rates. Some of this 
is obviously due to the beamed nature of the 
(variable) quasars, enhancing the variability both 
by boosting their brightnesses and shortening the 
variability timescales^ . Nonetheless, especially to- 
ward lower radio fiux density limits, a statistically 
significant study of variability is possible provided 
one starts with large enough samples. Large scale 
optical surveys like the SDSS provide the crucial 
radio source optical "environmental" information. 
As more SDSS data become available on the zero- 
dec strip, we will investigate this further. 
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^Even though blazars represent a very small fraction of 
AGN, they are much more common in variable samples. 
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. 10. — FZDA images of the extended subset of our variable sample (28 sources in total). The boxes 
x72" in size. The variable component is indicated by the green cross. 
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Table 1 

FZDA / FZDB sample statistics 



Sample 


Total Number 


Point Sources Only 


FZDA 


7343 


5550 


FZDB 


9086 


7286 


Combined 


6605 


5172 


FZDB not in A 


2481 


1846 



17 



Table 2 

FIRST Zero-Dec Variable Sample 



*J 1 11 li Cl 1 1 


RA (J2000) Dec (J2000) 
(core component) 


.^2002 

fm Tvl 


^ 1995 
fm Tvl 


FR^ 






J* 


C 


IDC 


RM*^ 


JOOOlll-002011 


uu 


ni 

Ui 


11 1 o 

ii.iy 


nn 

— UU 


on 
zu 


1 1 

ii. 


,D 


33.57 


23.28 


1 

i. 


,4D 


1 p; 

10 


.D 


0.518 


18, 


,56 


Q 


CL 


J000257-002447 


nn 
UU 


no 
Uz 


^7 1 a 
0( .iy 


nn 
— UU 


O/l 
Z4 


A 7 
4( , 


A 

,4 


112.93 


124.90 


n 
U, 


no 
.yz 


A 

4, 


o 
.z 




19, 


.52 


G 


CH 


J000549+005048 


UU 


Uo 


4y.yz 


nn 

UU 


ou 


48. 


1 

.i 


8.87 


5.13 


1 

i. 


77 


1 n 

iU 


1 

.i 




21, 


.59 


G 


PS 


J001158-000208 


nn 
UU 


1 1 

ii 


Q1 

Oo.oi 


nn 

— UU 


no 

uz 


no 
Uo, 


o 

,z 


5.95 


27.34 


n 

U, 


OO 

,zz 


/I Q 
4o 


o 
.z 




20, 


,98 


Q 


PS 


J001507-000801 


UU 


1 K 

io 


nv m 


nn 
— UU 


nc 
Uo 


ni 
Ui. 


Q 

,0 


10.54 


12.67 


n 
U, 


.00 


4, 


.6 


1.704 


18, 


.23 


Q 


PS 


J001611-001512 


nn 
UU 


1 a 
io 


1 1 n'7 
ii.U/ 


nn 
— UU 


1 c 
iO 


1 o 

iz. 


Q 

• O 


827.61 


1050.26 


n 

U, 


Qn 
.oU 


ii 


o 
.z 


1.575<^ 


19, 


.70 


Q 


PS 


J002354+001919 


nn 
UU 


zo 


04. OO 


nn 
UU 


1 n 

iy 


1 n 

iy. 


c 
,0 


4.75 


7.46 


n 

U, 


,DD 







.z 




22, 


.69 


G 


C.J 


J002738+000627 


nn 
UU 


Z 1 


QQ O^ 
00. zO 


nn 
UU 


UD 


07 
Z ( , 


.D 


4.56 


6.91 


n 
U, 


.Do 


0, 


A 

.4 




22, 


.02 


Q 


PS 


J003007-000007 


nn 
UU 


Qn 


n'7 on 

u ( .yu 


nn 
— UU 


nn 
UU 


n7 
U/ , 


.4 


50.52 


39.32 


1 

i. 


Qn 
.oU 


1 o 
iz 


1 

.i 




19, 


,55 


Q 


PS 


J003127+003959 


nn 
UU 


Q 1 

oi 


o'7 no 
z / .yz 


nn 
UU 


Qn 

oy 


p;n 

oy. 


,4 


9.58 


12.31 


n 

U, 


7n 


D 


.6 




21, 


,80 


Q 


PS 


J003536-000627 


UU 


60 


Qn 
oo.oU 


nn 
— UU 


UD 


07 
Z ( , 


o 
.z 


17.28 


14.52 


i. 


oi 
.zi 


D, 


• i 




18, 


.94 


G 


PS 


J003540-002529 


nn 
UU 


OO 


An OO 
4U. zz 


nn 
— UU 


op; 
zO 


OO 

zy, 


o 

,z 


4.01 


6.57 


n 

U, 


^?Q 

.Do 


D, 


.4 




21, 


.82 


G 


PS 


J003703+003537 


nn 

UU 


O ( 


nQ 

Uo.oO 


nn 

UU 


Qf^ 
OO 


Q7 
O ( , 


,4 


12.84 


15.28 


n 

U, 


,0D 


A 

4 


.0 




21, 


.51 


Q 


PS 


J004332+002459 


nn 
UU 




QO '71 
OZ. ( i 


nn 
UU 


O/l 
z4 


oy. 


o 
.0 


120.94 


108.53 


i. 


1 Q 

.io 


D, 


• i 


1.127'= 


19, 


.20 


Q 


PS 


J004819+001457 


nn 
UU 


A Q 
48 


1 O 1 1 

iy.ii 


nn 
UU 


1 /( 
i4 


P;7 
0/ , 


1 

.i 


98.76 


89.40 


1 

i. 


1 o 

.iz 


p: 
0, 


.0 


1.536<^ 


19, 


.92 


Q 


PS 


J005205+003538 


nn 
UU 


r^o 


UO.Oo 


nn 
UU 


Q 

oO 


QQ 

00, 


o 

,z 


34.58 


81.46 


n 

U, 


A Q 
,4o 


QQ 
08 


1 

. i 


0.399 


16, 


,09 


Q 


CJ 


J005212+000945 


nn 
UU 


OZ 


1 O /I '7 

iz.4 ( 


nn 
UU 


nn 

uy 


/I 

40, 


o 
.z 


11.92 


9.77 


i. 


O/l 
.z4 


0, 


.0 










PS 


J005225+002627 


nn 
UU 


OZ 


OK 

ZO.OI 


nn 
UU 


Ofi 
ZD 


OQ 

z8. 


n 
.U 


12.13 


9.98 


i. 


O/) 

.z4 


r 

0, 


o 
.z 




20, 


.75 


G 


PS 


J005717-002433 


nn 

UU 


( 


1 7 m 

i ( .Ui 


nn 

— UU 


O/l 

Z4 


QQ 
00, 




,z 


89.47 


114.74 


n 

u. 


70 

, 1 y 


1 1 

11 


.D 


2.790« 


19, 


.19 


Q 


CJ 


J010525+001121 


ni 
Ui 


UO 


zO.Oz 


nn 
UU 


1 1 
ii 


OI 
Zi, 


7 


33.97 


41.73 


n 
U, 


CQ 
.OO 


Q 

8, 


Q 

.6 




21, 


.45 


G 


PS 


J010745+003952 


ni 
Ui 


UI 


40. zi 


nn 
UU 


Qn 

oy 


p;o 
Oz, 


Q 

.8 


10.50 


8.88 


1 

i. 


OI 

.zi 


4, 


p; 
.0 










PS 


J010838+002814 


ni 
Ui 


no 
Uo 


OO.OD 


nn 
UU 


OO 
ZO 


1 A 

i4. 


,D 


5.33 


3.91 


1 

i. 


A n 

,4U 


4 


o 

.z 










PS 


J011106+000846 


ni 
Ui 


ii 


r\fy '7n 

ud. ( y 


nn 
UU 


no 
Uo 


A 

4D, 


.i 


5.25 


3.43 


i. 


C^7 

.0 / 


4, 


.y 




22, 


.70 


Q 


PS 


J011515+001248 


ni 
Ui 


1 

iO 


1 K 7Q 

iO. (o 


nn 
UU 


1 o 
iz 


48, 


p: 
.0 


46.84 


43.13 


1 

i. 


1 n 
.iU 


4, 


.4 


0.045 


14, 


.39 


G 


PS 


J012213-001801 


ni 
Ui 


OO 

zz 


1 Q no 

io.yz 


nn 
— UU 


1 

io 


ni 

Ui, 


n 

,U 


386.60 


331.60 


1 

i. 


1 

,io 


Q 

8 


Q 

.0 




20, 


,23 


Q 


CL 


J012528-000555 


ni 
Ui 


o 

zO 


oc Q 
zo.oO 


nn 
— UU 


UO 


00, 


Q 

.0 


1333.76 


1481.35 


n 
U, 


OI 

.yi 


4, 


Q 

.8 


1.076 


16, 


.47 


Q 


PS 


J012753+002516* 


ni 
Ui 


o'? 
z/ 


CO 7n 
Oo. /U 


nn 
UU 


op; 
zO 


1 

io. 


p: 
.0 


131.32 


90.08 


1 

i. 


/(C 
.48 


1 o 

iy 


1 

.i 




20, 


.76 


Q 


PS 


J013457+003942 


ni 
Ui 




/ .4z 


nn 
UU 


Qn 

oy 


A Q 

4o, 


n 
,U 


6.09 


2.87 


o 

z. 


1 Q 

,io 


Q 

8 


'7 
. ( 




22, 


,04 


Q 


PS 


J013815+001446 


ni 
Ui 


OO 


1 no 
iO.Uz 


nn 
UU 


1 A 

i4 


4D, 


r 
,0 


42.98 


51.76 


n 
U, 


C /t 

.84 


Q 

8, 


• i 




21, 


.64 


G 


CJ 


J015329-002214 


01 


53 


29.75 


—00 


22 


14, 


.3 


17.46 


14.98 


1, 


.19 


5, 


.8 




19, 


.07 


Q 


PS 


J015528+001204 


01 


55 


28.47 


00 


12 


04, 


,6 


17.23 


19.82 


0, 


,88 


4, 


.8 




22, 


,08 


G 


PS 


J015950-002407 


01 


59 


50.09 


-00 


24 


07, 


,2 


10.68 


12.63 


0, 


,86 


4, 


.3 




21, 


,50 


G 


PS 


J020141+003825 


02 


01 


41.04 


00 


38 


25, 


.5 


4.50 


6.35 


0, 


.73 


4, 


.7 










PS 


J020214-001748 


02 


02 


14.30 


-00 


17 


48, 


,3 


75.12 


60.00 


1, 


,27 


11 


.5 




21, 


,36 


G 


PS 


J020234+000301 


02 


02 


34.32 


00 


03 


01, 


,7 


30.54 


39.41 


0, 


,79 


10 


.5 




18, 


.42 


G 


CL 


J020928-001224 


02 


09 


28.85 


-00 


12 


24, 


,9 


4.10 


2.51 


1, 


,69 


4, 


.2 


0.152 


16, 


.02 


G 


PS 


J021202-002750 


02 


12 


02.13 


-00 


27 


50, 


,1 


45.14 


52.84 


0, 


,87 


6, 


.6 










CL 


J021301-001815 


02 


13 


01.13 


-00 


18 


15, 


,0 


41.55 


48.44 


0, 


,87 


6, 


.5 




21, 


,98 


G 


PS 


J021553+001826 


02 


15 


53.65 


00 


18 


26, 


,9 


30.90 


35.82 


0, 


,88 


5, 


.6 




19, 


,51 


G 


PS 


J021755-000935 


02 


17 


55.99 


-00 


09 


35, 


,9 


3.96 


5.85 


0, 


,70 


4, 


.4 




19, 


,86 


G 


CJ 


J021840-001515 


02 


18 


40.55 


-00 


15 


15, 


,9 


10.06 


13.17 


0, 


,78 


6, 


.4 


1.171 


18, 


,81 


Q 


cx 


J022624+000746 


02 


26 


24.61 


00 


07 


46, 


,2 


43.84 


48.74 


0, 


,91 


4, 


.4 




22, 


,56 


G 


PS 


J023105+000843 


02 


31 


05.59 


00 


08 


43, 


,5 


47.60 


53.86 


0, 


,90 


5, 


.0 


1.338 


20, 


,00 


Q 


PS 


J025321+000559 


02 


53 


21.04 


00 


05 


59, 


,9 


151.89 


97.64 


1, 


.57 


22, 


.0 




22, 


.26 


G 


PS 



Table 2 — Continued 



Source / 


RA (J2000) Dec (J2000) F2002 


-fl995 FR*^ 




1'= m'= RM^ 


Component 


(core component) [mJy] 


[mJy] 







J025333+002431 


02 


53 


33 


.66 


00 


24 


31, 


.7 


16, 


.86 


14 


.86 


1 


.15 


4. 


,6 








PS 


J025404-002628 


02 


54 


04, 


.60 


-00 


26 


28, 


,7 


4, 


,39 


2, 


.71 


1 


.67 


5, 


,2 




21.97 


G 


PS 


J025859+003618 


02 


58 


59, 


.65 


00 


36 


18, 


.3 


29, 


.45 


23 


.75 


1 


.26 


9. 


,4 




22.28 


Q 


PS 


J025928-002000f 


02 


59 


28, 


.52 


-00 


20 


00, 


,1 


243, 


.47 


225 


.85 


1 


.09 


4. 


,3 


2.001 


17.34 


Q 


PS 


J030137+001746 


03 


01 


37, 


.58 


00 


17 


46, 


,1 


6, 


,86 


4, 


.45 


1 


.58 


7, 


,0 




21.47 


G 


CJ 


J030317-002644 


03 


03 


17, 


.01 


-00 


26 


45, 


,0 


7. 


,41 


3, 


.68 


2 


.06 


4. 


,9 








CJ 


J030702+000651 


03 


07 


02, 


.04 


00 


06 


51, 


.9 


5. 


,84 


8, 


.08 





.74 


5. 


,4 








PS 


J030834+003303 


03 


08 


34, 


.31 


00 


33 


03, 


,7 


4, 


,16 


1, 


.83 


2 


.35 


6, 


,4 


0.031 


14.91 


G 


PS 


J030933-001901 


03 


09 


33, 


.27 


-00 


19 


01, 


,2 


6. 


,31 


4, 


.58 


1 


.41 


4. 


,6 








PS 


J031006+001549 


03 


10 


06, 


.69 


00 


15 


49, 


.9 


1. 


,36 


5, 


.55 





.26 


11, 


.1 


0.109 


18.54 


G 


PS 


J031118+000816 


03 


11 


18, 


.51 


00 


08 


16, 


.4 


29, 


.30 


36 


.54 





.81 


9, 


,1 








PS 


J031202-000442 


03 


12 


02, 


.50 


-00 


04 


42, 


,5 


6. 


,42 


8, 


.59 





.77 


5. 


,2 


0.038 


13.50 


G 


PS 


J031345-000720 


03 


13 


45, 


.04 


-00 


07 


20, 


.3 


3. 


,63 


6, 


.06 





.62 


6. 


,1 


2.519 


20.04 


Q 


PS 


J031353-000004 


03 


13 


53, 


.47 


-00 


00 


04, 


,1 


10, 


,13 


7, 


.03 


1 


.47 


7, 


,7 




21.16 


G 


PS 


J031357+003506 


03 


13 


57, 


.10 


00 


35 


06, 


.9 


16, 


.54 


20 


.23 





.83 


6. 


,7 








PS 


J031452+001346 


03 


14 


52, 


.08 


00 


13 


46, 


.5 


6. 


,77 


10 


.05 





.69 


7. 


,3 


3.202 


19.54 


Q 


PS 


J031609+000107 


03 


16 


09 


.55 


00 


01 


07, 


,9 


4, 


,02 


6, 


.03 





.69 


4, 


,6 








PS 


J031634-002039 


03 


16 


34, 


.96 


-00 


20 


39, 


,6 


18, 


.17 


15 


.63 


1 


.18 


5. 


,5 








PS 


J031814-002948f 


03 


18 


14, 


.43 


-00 


29 


48, 


.9 


113, 


.20 


93 


.91 


1 


.22 


9. 


,7 




21.77 


G 


PS 


J032007+000354 


03 


20 


07, 


.01 


00 


03 


54, 


,0 


12, 


,17 


10 


.20 


1 


.22 


5, 


,0 




22.10 


G 


PS 


J212000-001159 


21 


20 


00, 


.72 


-00 


11 


59, 


,5 


42, 


.94 


48, 


.33 





.90 


4. 


,9 








PS 


J212447+000747 


21 


24 


47, 


.34 


00 


07 


47, 


.6 


5. 


,79 


3 


.82 


1 


.56 


5. 


,3 








PS 


J213206+003520 


21 


32 


06 


.15 


00 


35 


20, 


.1 


13, 


,65 


11 


.95 


1 


.16 


4, 


,1 








PS 


J213638+004154f 


21 


36 


38, 


.57 


00 


41 


54, 


.3 


4136, 


.68 


3546 


.71 


1 


.18 


8. 


,3 


1.932« 


16.79 


q 


PS 


J213748+001219 


21 


37 


48, 


.43 


00 


12 


19, 


.9 


41, 


.69 


36 


.02 


1 


.17 


7. 


,0 


1.666<^ 


17.92 


q 


PS 


J213804+000714 


21 


38 


04, 


.06 


00 


07 


14, 


,7 


10, 


,31 


16 


.33 





.64 


11, 


.6 








CJ 


J214138+000319 


21 


41 


38, 


.55 


00 


03 


19, 


.8 


11, 


.08 


14 


.37 





.79 


6. 


,8 








PS 


J214324+003502 


21 


43 


24, 


.37 


00 


35 


02, 


.7 


36, 


.03 


45, 


.22 





.81 


9. 


,7 


2.030^ 


19.37 


q 


PS 


J214419+002055 


21 


44 


19, 


.88 


00 


20 


55, 


,8 


10, 


,34 


7, 


.31 


1 


.44 


7, 


,5 








PS 


J214613+000930 


21 


46 


13, 


.31 


00 


09 


30, 


.8 


9. 


,94 


7, 


.88 


1 


.29 


5. 


,6 








PS 


J214807-000630 


21 


48 


07, 


.73 


-00 


06 


30, 


.6 


68, 


.57 


90 


.18 





.77 


12, 


.7 








CJ 


J214811-001230 


21 


48 


11, 


.48 


-00 


12 


30, 


.6 


7, 


,76 


9, 


.95 





.80 


4, 


,7 








CJ 


J215349+003119 


21 


53 


49, 


.75 


00 


31 


19, 


.6 


82, 


.50 


103 


.65 





.81 


10, 


.4 








PS 


J215353-001339 


21 


53 


53, 


.89 


-00 


13 


39, 


.5 


13, 


.48 


17 


.27 





.79 


7. 


,3 








PS 


J215359+004412 


21 


53 


59, 


.83 


00 


44 


12, 


,6 


9, 


,76 


6, 


.73 


1 


.48 


7, 


,8 


1.030« 


19.07 


q 


PS 


J215733-000340f 


21 


57 


33, 


.66 


-00 


03 


40, 


.5 


5. 


,36 


3, 


.60 


1 


.53 


4. 


,9 








PS 


J215949+005146 


21 


59 


49, 


.91 


00 


51 


46, 


.7 


12, 


.29 


9, 


.21 


1 


.36 


7. 


,8 








PS 


J215954-002150 


21 


59 


54, 


.46 


-00 


21 


50, 


.1 


4, 


,07 


2, 


.46 


1 


.71 


4, 


,7 


1.960« 


16.98 


q 


PS 


J220017-000133 


22 


00 


17, 


.37 


-00 


01 


33, 


.6 


7. 


,84 


6, 


.15 


1 


.30 


4. 


,9 








CL 


J220755-000215 


22 


07 


55, 


.25 


-00 


02 


15, 


.0 


78, 


.40 


61 


.93 


1 


.28 


11, 


.9 








CJ 


J220822+002352 


22 


08 


22 


.88 


00 


23 


52, 


,7 


4, 


,62 


2, 


.42 


1 


.97 


6, 


,1 








PS 


J221001-001309f 


22 


10 


01, 


.82 


-00 


13 


09, 


.9 


125, 


.83 


115 


.31 


1 


.10 


4. 


,8 








PS 


J221031-001356 


22 


10 


31, 


.46 


-00 


13 


56, 


.1 


14, 


.06 


12 


.00 


1 


.19 


5. 


,2 








PS 


J221909+003113 


22 


19 


09, 


.40 


00 


31 


13, 


.4 


9. 


,13 


11 


.70 





.80 


5. 


,3 








CJ 



19 



Table 2 — Continued 



Source / 


RA (J2000) 


Dec (J2000) 


-f2002 


-^^1995 


FR*^ 


S'' 








RM'^ 


Component 


(core component) 


[mJy] 


[mJy] 
















22 20 36.33 


00 33 34.2 


1 97 


1 9 7fi 
iz. r D 


1 99 
l.ZZ 


U.O 








PC! 

1^0 


79991'?^ nminn 

jzzzj_oo UUllUU 


22 21 35.00 


-00 11 00.1 


14 11 


8 R4 

O.U4 


1 .uu 


1 9 fi 
1 z .u 










Tf)f)f)f)or I nrii roc 


22 22 35.87 


00 15 36.7 


4fi n4 


1^1 1 


n 01 

U.y 1 


4 

^.u 








P7 


T9997n44-nfl4'i1 7 


22 27 04.24 


00 45 17.5 


8 13 

O. lO 


■i "^7 


1 4Q 
1 .iy 


7 n 

1 .u 








PS 


T99979fi-l-nm (^'^Q 


22 27 26.53 


00 10 59.3 


fi n^ 


4 ^i9 


1 37 

1 .0 ( 


4 9 

4:. Z 








PS 


T99979Q-l-flflfl'i99 


22 27 29.08 


00 05 22.2 


77 84 

i i .OH; 


Q1 fi4 
y 1 


n 8fi 


7 4 
1 .H 


1 'iin'^ 

1 .U lU 


18 8fi 
10. ou 


Q 


r!T, 


7999744-1- flfl'^d'ifl 


22 27 44.58 


00 34 50.6 


99 1 3 


3n 113 


n 74 


1 9 1 

IZ. 1 


1 .U^U 


1Q 14 
ly . 1^ 


Q 


PS 


T9997^8+nn'?7n'=l 


22 27 58.13 


00 37 05.2 


fi8 RQ 


QQ 1 9 
y y . 1 z, 


n 70 


1 7 3 
1 1 .o 








PS 


T99'^n47-l-nn97iifi 


22 30 47.46 


00 27 56.3 


S ii4 


1 89 
lo.oz 


n (v\ 


114 
11.^ 








1 


T994994-l-fin'i'i1 9 


22 42 24.14 


00 55 13.0 


1 7 71 


1 3 44 


1 34 
1 .o^ 


y .u 










T994'^'^1 —Dm ')X\ 


22 43 31.94 


-00 12 33.1 


1 3 8^ 

lO . 


91 40 


n fifi 

u .uu 


1 4 8 

14:. O 


Z .U^U 


18 411 
10. ^u 


Q 


PS 




22 44 48.11 


-00 06 19.8 


8 94 




1 l=.1 

1 . U 1 


7 3 








PS 


T994fi97— nni 91 4^ 


22 46 27.68 


-00 12 14.2 


8"! 'I'l 


•ifi nn 


1 .uu 


9n Q 

zu. y 








PS 


T9947'^n+nnnnnfif 


22 47 30.19 


00 00 06.1 


4fi4 ^9 

4:U4: . (JZ; 


1 83 71 

1 OO . ( 1 


z . uu 


43 8 

40.0 


n nQ4<* 

u.uy^ 


1 8 i^n 

10. uu 


Q 


PS 


T99/lCl99-J-nm Sfl/l 


22 49 22.28 


00 18 04.6 


1 n 41 

lU.'ll 


S 1^1 


1 9^ 
1 .zu 


U.O 








1 


T99iiSii9— nni S'=i7 


22 58 52.94 


-00 18 57.3 


4^^ 


1 98 
1 .zo 


9 sn 

z.ou 


fi n 

u.u 








p*^ 




23 01 57.85 


00 03 52.0 


3 7^1 


7 n4 


n 

u . uu 


7 8 

( . O 




93 93 




p*^ 


T9'^n'^i 4-i-nnnnti9 


23 03 14.85 


00 00 52.2 


% 79 


1 79 
1 . J z 


9 94 

Z.Z1 


8 
U.O 








p*^ 


T9'^nfi'i'i-i-nn'^fi'^8 


23 06 55.16 


00 36 38.1 


1 9 QQ 
1^ . 


1 1 Q 
i(j. 1 y 


n 87 


4 3 








p^ 


T9'?1 ^41 -t-nn9Q'?R 


23 15 41.67 


00 29 36.6 


91 Q"? 


1 Q4 


1 fin 

1 .uu 


1 '^ 7 

lu. ( 




9n 81 

ZU.Ol 


n 


pcl 


79*^1 i^t^s nni9n'^ 


23 15 58.64 


-00 12 05.5 


% 1^9 


(\ 97 
O.Z / 


n t^s 

u.uo 


7 n 




99 8'^ 
zz. ou 


p 

V I 


1 kJ 


T9'?i SAii— nnn7'=i4 


23 18 45.81 


-00 07 54.7 


O.VJO 




n 

U.Uy 


U.O 


n 8fi7 
u.ou / 


1 Q 44 
ly .44 




p*^ 




23 18 56.66 


00 14 37.7 


90 1 (\ 

ZU. lU 


1 8 '^4 


1 1 9 

1 . 1 z 


4 9 
4. z 


n n'^n 

u.uou 


1 9 7fi 
IZ. ( u 


c 


pc 


T23iqin-i-nni85q 


23 19 10.33 


00 18 59.0 


28.90 


33.36 


0.88 


5.6 




22.60 





PS 


7232038-1-003139 


23 20 38.01 


00 31 39.8 


72.69 


82.95 


0.89 


5.7 


1.911 


19.03 




PS 


7232323+003327 


23 23 23.95 


00 33 27.5 


16.73 


12.52 


1.36 


9.2 








PS 


T233448— 001400 


23 34 48.06 


-00 14 01.0 


24.85 


29.64 


0.85 


6.8 




21.56 





PS 


7933829-1-001 14fi 


23 38 22.35 


00 11 46.6 


18 1Q 

J.O. J- 1/ 


15.44 


1.20 


fi 




22.24 





PS 


79338^9+004843 


23 38 52.46 


00 48 43.5 


11.32 


Q 1 8 


1.26 


5 Q 








PS 


J234623+004301 


23 46 23.73 


00 43 01.1 


13.49 


11.79 


1.17 


4.6 


2.861 


19.09 


Q 


PS 


J234939-001315 


23 49 39.90 


-00 13 15.2 


8.46 


6.97 


1.24 


4.2 


1.267 


20.22 


Q 


PS 


J235022+001232 


23 50 22.41 


00 12 32.4 


3.92 


2.03 


2.00 


5.1 








PS 


J235050-002848 


23 50 50.72 


-00 28 48.7 


13.12 


11.57 


1.15 


4.1 




19.62 


G 


CL 


J235409-001948f 


23 54 09.18 


-00 19 48.1 


384.58 


344.79 


1.13 


6.1 


0.462 


17.93 


Q 


PS 


J235823+000213 


23 58 23.91 


00 02 13.2 


8.03 


10.67 


0.77 


5.7 








PS 


jooisoo+oooaia*^ 


00 18 00.79 


00 03 17.9 


20.76 


18.82 


1.12 


4.2 








HS 


J003246-001917^ 


00 32 46.02 


-00 19 17.8 


81.29 


60.80 


1.35 


14.5 








HS 


J212058+000612h 


21 20 59.00 


00 06 12.7 


31.17 


28.05 


1.13 


5.3 








HS 


J212955+000758'^ 


21 29 55.68 


00 07 59.0 


1.17 


3.03 


0.42 


4.5 








CL 


J235828+0039341' 


23 58 28.77 


00 39 34.1 


2.48 


4.16 


0.62 


4.3 








CL 



"Corrected flux density ratio, sec Eqn. a. 
'^Significance (in units of a rms) as defined by: S 



x(FR+i) 

SDSS redshift, r-band magnitude, and optical identification. 



FR- 



, using the % definition from Eqn. 4. 

except where noted. Sources in the hst with RA 
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between 21^ and 23*^ are not covered by SDSS DRl.. , , , ^„ 

Radio morphology: PS=pomtsource, UJ=core-jet, CL=core + lobes, CA=coniplcx, CH=core + halo, and 

''^^R69i&\'Fl!'tfefiiSPft|jm htcratui-c. If the ID is in lower case then the optical data arc from the literature as well. 
^NUy^^l}^ the FZDB-NVSS variability criterion. Only sources brighter in the FZDB epoch were considered. 



h 



Variability of this source is spurious (see Sect. 4) 
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Table 3 

apm match statistics, 3" search 

RADIUS 



Sample** 


Detected 


Rate 


V CXI. lex UlC 


66 


53.7% 


Control 1 


17 


13.8% 


Control 2 


29 


23.6% 


Control 3 


21 


17.1% 


Control 4 


21 


17.1% 


Control 5 


20 


16.3% 


Control 6 


27 


22.0% 


Control 7 


23 


18.7% 


Control 8 


17 


13.8% 


Control 9 


21 


17.1% 


Control 10 


22 


17.9% 


Control 11 


23 


18.7% 


Control 12 


21 


17.1% 


Mean Control 


21.8±3.5 


17.7±2.8% 



"^AU samples contain 123 sources 
Table 4 

sdss match statistics, 3" search radius 



Sample** 


Detected 


Rate 


Quasars'' 


Galaxies^ 


Variable 


64 


82.0% 


53.1% 


46.9% 


Control 1 


28 


37.3% 


17.9% 


82.1% 


Control 2 


30 


40.0% 


13.3% 


86.7% 


Control 3 


21 


28.0% 


14.3% 


85.7% 


Control 4 


36 


48.0% 


13.9% 


86.1% 


Control 5 


26 


34.7% 


26.9% 


73.1% 


Control 6 


32 


42.7% 


18.8% 


81.2% 


Control 7 


29 


38.7% 


24.1% 


75.9% 


Control 8 


33 


44.0% 


9.1% 


90.9% 


Control 9 


23 


30.7% 


21.7% 


78.3% 


Control 10 


30 


40.0% 


16.7% 


83.3% 


Control 11 


31 


41.3% 


29.0% 


71.0% 


Control 12 


22 


29.3% 


18.2% 


81.8% 


Mean Control 


28.4±4.6 


37.9±6.2% 


18.7±5.9% 


81.3±5.9% 



'^Only 75 out of 123 variable sources are covered by the SDSS DRl 
data. The rates have been corrected accordingly. 

''Relative fraction of quasars and galaxies, based on the SDSS mor- 
phology. 
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Table 5 

Mean Optical colors of Variable and 
NoN- Variable Samples 



Sample 


(U-R) 


(R-Z) 


Quasars - variable 


0.89±0.13 


0.40±0.07 


non- variable 


1.37±0.17 


0.60±0.G8 


nominal offset^ 


-0.48 


-0.20 


significance'' 


2.3 


1.8 


GalcLxies - variable 


1.81±0.21 


1.09±0.11 


non- variable 


2.28±0.09 


1.20±0.05 


nominal offset*^ 


-0.47 


-0.11 


significance'' 


2.0 


0.9 



'^Between non- variable and variable samples. 

''Defined as nominal offset divided by the rms 
of the uncertainties. 
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